Entangling remote nuclear spins linked by a chromophore 
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Molecular nanostructures may constitute the fabric of future quantum technologies, if their degrees of 
freedom can be fully harnessed. Ideally one might use nuclear spins as low-decoherence qubits and optical 
excitations for fast controllable interactions. Here, we present a method for entangling two nuclear spins 
through their mutual coupling to a transient optically-excited electron spin, and investigate its feasibility 
through density functional theory and experiments on a test molecule. From our calculations we iden- 
tify the specific molecular properties that permit high entangling power gates under simple optical and 
microwave pulses; synthesis of such molecules is possible with established techniques. 



Molecules are promising building blocks for quantum 
technologies, due to their reproducible nature and ability to 
self-assemble into complex structures. However, the need 
to control interactions between adjacent qubits represents 
a key challenge EH?). We here describe a method for op- 
tical control of a nuclear spin-spin interaction that presents 
several advantages over conventional NMR quantum pro- 
cessors: the spin-spin interaction can be switched, the gates 
are faster, and the larger energy transitions facilitate polar- 
ization transfer onto the nuclear spins. After explaining 
the theory behind our method, we present an experimental 
study of a test molecule, and show with density functional 
theory that an entangling gate could be achieved. 

We consider two nuclear spin qubits labelled n and n' 
and one mediator e. The mediator possesses a paramag- 
netic excited state \e) with spin one character and a dia- 
magnetic, spinless ground state |0). The two nuclear qubits 
do not interact with each other directly, but both couple to 
the excitation via an isotropic hyperfine (HF) coupling with 
generally unequal strengths A and A', see Fig. [I] (a). The 
Hamiltonian in a magnetic field is given by (h = 1): 

H = -(O n S Zin - (D n iS z , n > + \e) ((O e S Zie + (Do) (e\ 

+ \e)(AS n 'S e +A f S^'S e +DSl e )(e\ , (1) 

where S ZjI - and S; are the Pauli spin operators and Cty de- 
notes the respective Zeeman splittings (i = n, n',e). D is the 
zero-field-splitting (ZFS). (Do denotes the optical frequency 
corresponding to the creation energy of the excitation. 

Let us first analyze the case of a symmetric (homonu- 
clear) system with (D n = (D n > and A = A 7 . Using degener- 
ate perturbation theory and assuming that the electronic 
Zeeman splitting is much larger than that of the nuclei, 
(D n <C (D e and the ZFS, D <C (D e , we obtain an effective 
Hamiltonian by approximating H sym := (e\H\e) in the fol- 



lowing way: 



H sym « V^diag^i , . . . ,Ei 2 )V =: # sy m,eff 



(2) 

where V = ( | ) | Xfo) • • • | y/12) ) is the matrix of the approxi- 
mate eigenvectors up to first order and the E\ (i — 1 , . . . , 12) 
are the eigenenergies up to second order. This reveals that 
the time evolution of the entire system can be decomposed 
into three different unitary actions on the nuclear spin, one 
for each spin projection of the excitation |7}) (i = — ,0,+): 



U(t) w exp(-i# sym ,eff0 = U-(t)®U (t)®U+(t) . (3) 

Moreover, only the states |7] ]/[) and |7] are coupled, 
with a matrix element af. 
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This coupling can take an initial product state to an en- 
tangled state of nuclear spins and we can quantify the de- 
gree of entanglement by employing the entangling power 
ei of a unitary U(. e\ is the mean linear entropy produced 
by the unitary acting on a uniform distribution of all pure 
product states. A maximally entangling gate has e\ = | |[5l. 
For Eq. ([3]) we get 

ei(t) = i(3 + cos(2fl£f)) sin 2 (ait) for / = - ,0,+ . (5) 

Hence U(t) is maximally entangling at odd integer multi- 
ples oft—^-. eo(t) w for t < ^- since ao <C a± (see 

Fig.0. 

A protocol for creating and preserving entanglement of 
the nuclear qubits is thus to start in the singlet ground state 
of the electron and to initialize the nuclear spins, e.g. in 
the state 1 Polarizing nuclear spins in an NMR exper- 
iment is very challenging, however in this case one could 
initialize the nuclear spin by transferring the electronic spin 
polarization. A laser pulse then creates the electron triplet. 
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FIG. 1 : (Color online) (a) A triplet state of an electron and hole 
couples two nuclear spin qubits. The triplet can be |7+) = | tA)» 
\To) = j^{\ ieW + I teU)) or |r_) = | ieiu). (b) An example 
molecular embodiment of such a model system is bis-diethyl mal- 
onate fullerene. (c) Entangling power of the time evolution oper- 
ators U+ (solid curve) and Uq (dashed curve); here a+ « 32^o. (d) 
Maximal attainable entangling power ra_ for U- as the symmetry 
of the two nuclear spins is reduced. Ai = fe^^J ; A2 = ; 
A = 2.5 MHz; D = -296 MHz; co e = 9.6 GHz;"co„ = 3.7 MHz. 

Our experiments (see below) revealed that the state of the 
excitation depends on the orientation of the molecule to 
the static magnetic field. As a start we assume it is in |7b). 
With a strong microwave pulse we flip the state of the exci- 
tation, say to the state, and then wait a time t — ^-.A 
second microwave pulse is then used to flip the excitation 
back into |7o). If the lifetime of the excitation T is bounded 
by 

(6) 

the nuclear spins remain entangled for the duration of nu- 
clear spin coherence (26), because the electronic ground 
state is spin silent. 

Such a scheme could be demonstrated using an ensem- 
ble of molecules, so long as we address two challenges. 
First, spontaneous optical decay must be taken into ac- 
count. In general, the decay process gives rise to distin- 
guishable photons, destroying coherence between the var- 
ious spin states. However, in the limit of Eq. ([6]), the 
states I To ti) and \Tq |t) decay producing indistinguish- 
able photons, and, therefore, nuclear spin coherence in this 
space will survive the excitation decay. Specifically, ao is 
much smaller than the natural linewidth of the optical tran- 
sition (unless A assumes a very large value). Second the 



qubits must be individually addressable to verify entan- 
glement generation with state tomography. This requires 
an asymmetry in (O n i . Generally, the eigenspectrum of an 
asymmetric molecule is not degenerate, and this tarnishes 
our scheme since e\ is then close to zero for all times. 
However, perturbation theory shows that under the condi- 
tion A' — A « ±2(a) n > — (On), we can recover a degenerate 
eigenspectrum for the \T±) states even for an asymmetric 
molecule. This leads to a very similar behaviour for the 
entangling power as in Eq. ([5]), see Fig.[T](d). 

There are many molecules possessing optically-excited 
triplet states which couple to nearby nuclear spins. These 
include pentacene, pyrazine, and 13 C6o ll6UIOl. We here 
determine the suitability of using a functionalized C60 
molecule for this scheme. Fullerenes possess a reason- 
able absorption coefficient at 532 nm and a highly effi- 
cient intersystem crossing (ISC) to a hyperpolarized triplet 
state HUE). Moreover, the number of spin active nuclei 
in these molecules is small, suppressing unwanted deco- 
herence fT3l . The HF couplings between 13 C atoms on the 
C60 cage and the photoexcited triplet are anisotropic with 
the largest HF term being 167.6 MHz (TUHEl. However, 
since the cage distortion in the excited state introduces a 
random element to the hyperfine coupling, we shall con- 
sider a functionalized fullerene whose molecular axis en- 
sures a well-defined interaction. 

Ultimately, an entangling operation would require two 
13 C spins, such as those of the trans- 1 bis-adduct [FIG. [I] 
(a,b)]. Individual addressibility could be achieved with dif- 
ferent chemical shifts of the two nuclei (which typically do 
not exceed 200 parts per million [ 13)) by making a small 
modification to one of the groups. In this case Fig. [T](d) 
reveals that with equal coupling constants A and A' an op- 
timal entangling operation is still achievable. Here we take 
the first step towards realizing this scheme by studying the 
mono-functionalized analogue to determine the key param- 
eters: hyper-polarization populations, triplet lifetime, and 
the HF tensor with a 13 C labeled methano-carbon of the 
diethyl malonate mono-adduct (DEMF) [FIG. [2] (inset)]. 

EPR measurements on the photo-excited paramag- 
netic state of spin labeled DEMF were performed at X- 
band (9.47 GHz) on a Bruker Elexsys580e spectrometer, 
equipped with a helium-flow cryostat. Photo-excitation 
was performed using a 532 nm pulsed Nd-YAG laser with 
10 Hz repetition rate; 7 ns pulse length; 5 — 10 mJ/pulse. 
The samples were prepared as 4.3-10 _4 M solutions in 
toluene-dg, were deoxygenated using five cycles of freeze- 
pump-thaw, flame sealed under vacuum and flash frozen in 
liquid nitrogen. 

Upon excitation of the fullerene cage to the singlet ex- 
cited state the system relaxes through ISC to the lowest 
triplet state which decays back to the singlet ground state in 
90—100 (as at 20 — 50 K, as measured by transient absorp- 
tion. The ISC relaxation mechanism preferentially popu- 
lates the triplet sublevels (T x ^ z ) according to the molecular 
symmetry. In the presence of a magnetic field the triplet 
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sublevels mix in an orientation dependent manner. This 
means in principle polarization could be achieved by ori- 
enting the molecules in a crystal. The electron spin echo- 
detected spectra revealed no change in the line shape as 
a function of excitation delay time or temperature in the 
range 5 — 50 K. This is in agreement with previous studies 
showing that at low temperatures, spin lattice relaxation ef- 
fects are negligible |[T6l . EPR line shape simulations were 
performed using the software package Easy spin |[T71l . The 
g tensor and ZFS parameters were consistent with pre- 
vious work [g xx = 2.0006,^ = 2.001 15, g zz = 2.00215, 
D = -296 MHz and E = -6 MHz] T\M . 

EPR line shape analysis yields only relative differ- 
ences in populations pi of the triplet sublevels (e.g. (p x — 
Py)/(Py ~ Pzj) 033 - To determine absolute populations 
and triplet decay rates for the three sub-levels we inves- 
tigated the relaxation kinetics, as shown in Fig. [2j In our 
analysis we neglect any decay between the sub-levels and 
consider separate decay rates k x ^ z from each sub-level to 
the singlet ground state. These decay rates and the initial 
populations n x ^ z are each properties of the triplet eigen- 
states at zero applied magnetic field (T x ^ z ), and so can 
be used to determine the resulting decay rates and pop- 
ulations at any applied magnetic field by writing the in- 
field eigenstates (r +; o,-) in the appropriate basis [20]. Us- 
ing this approach, we simultaneously fit the decay traces 
of flash delay (hv-T-n /2-T-;r-echo) and inversion recovery 
(/iv-t-;r-T-;r/2-T-;r-echo) experiments at four fields, using 
one set of parameters: p x : p y : p z = 0.46 : 0.54 : 0.00 
and x x = (0.50 ± 0.02) ms, x y = (0.58 ± 0.02) ms, x z = 
(0.020 ±0.003) ms (20 K), where x { = 1/jfc*. The transient 
populations and triplet decay rate were found to show lit- 
tle temperature dependence below 50 K. Importantly, the 
observed hyperpolarization would allow our electron me- 
diated gate to work well, and could be used to initialize the 
nuclear spins. 

The coupling of the triplet electron and the 13 C spin la- 
beled methano-carbon was measured by Mims electron- 
nuclear double resonance (ENDOR) EH at six field lines, 
see Fig. [2] Through simulations the hyperfine tensor was 
found to be axial, A xx = A yy = 2.6 ± 0.1, A zz = 0.74 ± 
0.02 MHz, though a 20% strain in A across the sample was 
required to fit the ENDOR linewidth at all fields (orienta- 
tions) ifTTll . This HF tensor is much smaller than that of 
certain 13 C atoms in the fullerene cage fTUl . 

In order to understand the weak HF coupling found 
from ENDOR measurements, and to predict the appropri- 
ate coupling of bis-DEMF, we turn to first-principles cal- 
culations with density-functional theory (DFT). For this 
purpose we considered optical excitations to spin triplet 
states in mono- and bis-diethyl malonate functionalized 
C6o molecules. The calculations were performed within 
the local density approximation to DFT [see supplemen- 
tary material for details]. To determine the ground-state 
structural configuration of the mono-diethyl malonate C6o 
we relaxed three different conformations with symmetries 
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FIG. 2: (Color online) (a) Decay traces and simulations of the 
flash delay (FD) and inversion recovery (IR) experiments at 20 K. 
Upper inset: DEMF molecule. Lower inset: electron spin echo 
field sweep of DEMF at 20 K. (b) Mims ENDOR simulations 
and experimental data acquired at 20 K at fields A-F from (a). 
The largest peak, around 3.9 MHz, was found to correlate to the 
13 C-7o interaction, while the weaker peak, at 2.2 MHz, related to 
2 H-7o interaction from the solvent (toluene-d8). Also observable 
in the ENDOR spectrum is the field dependent component of the 
hyperfine from the l3 C-T± triplet levels. 



C s , C2, or C2v The molecules with C2 or C s symmetry were 
found to be almost degenerate in energy (within 5 meV), 
and more stable than the molecule with C2 V symmetry by 
0.32 eV. The excitation to the spin-triplet state was stud- 
ied in detail for the molecule with C2 symmetry. By carry- 
ing out analogous structural relaxations for bis-diethyl mal- 
onate we found that the lowest energy configuration carries 
C2 symmetry. The spin-triplet excited states were calcu- 
lated by relaxing the electronic structure of each molecule 
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FIG. 3: (Color online) (a) Spin density in the lowest triplet ex- 
cited state of mono- and bis-diethyl malonate fullerenes. The 
chosen charge isosurface has an isovalue of 0.07 electron/A 3 and 
the overlaid color indicates the magnitude of the spin density. The 
spin density is mostly gathered around the fullerene cage. Images 
rendered using Gabedit (22). 

with the total spin constrained to S z = 1 E3l l24l . Fig. [3] 
shows that there is significantly more spin polarization on 
the cage than on the adducts. The calculated HF coupling 
constants of the 13 C nuclei in mono- and bis-dethyl adduct 
C6o are 4 MHz and 2.2 MHz, respectively. These values 
are in good agreement with our relatively small measured 
HF coupling. 

Our experiments on the DEMF molecule provide exam- 
ple triplet relaxation times and hyperpolarization of triplet 
sub-levels from which we can deduce the strength of HF 
coupling required to create significant nuclear entangle- 
ment [Fig. |4). When the laser excitation only produces 
partial electronic polarization, we need to adjust the pre- 
viously described protocol for creating and preserving en- 
tanglement. We consider the special situation of having all 
molecules oriented along the z-axis with parameters as ex- 
tracted from the experiments on DEMF, for which the pop- 
ulation would be almost equally distributed between \T+) 
and 1 71), and the optical decay from |7b) is much faster 
than from \T±). In this case, we can apply our protocol 
twice in succession since the entangling power is periodic 
in time: first for the population initially in \T+) and, af- 
ter waiting for |7b) to have emptied out, a second time for 
the population initially in |71). Given a two nuclear spin- 
system exhibiting a hyperfine coupling similar to that ob- 
served in our test molecule (3 MHz), one could perform 
an operation resulting in a state with entanglement of for- 
mation ^0.5. The hyperfine coupling for the bis-adduct 
fullerene could be enhanced by using an optimized func- 
tional group. In fact, the hyperfine coupling can be re- 
garded as a tunable parameter, and using the results of 
Fig. [4] an optimal value can be chosen to maximize the 
entangling power. In Fig. [4] the entanglement goes down 
for large coupling constants, because in this case the right- 
hand side of Eq. ([6]) becomes small enough to produce 
partially distinguishable photons in the decay process. 

Our analytic, numerical and experimental results have 
demonstrated the feasibility of manipulating nuclear spin 
qubits in molecular structures by harnessing optically ex- 
cited electron spin states. We identify the required molec- 
ular characteristics; it is feasible to synthesise suitable 
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FIG. 4: Entanglement of formation [25] of the nuclear spins after 
applying our scheme and when the excitation has decayed as a 
function of the HF strength for different initial polarization of the 
excitation. The nuclear spins for both curves are assumed to be 
initially in the state | J/f); parameters as in Fig.[l] For a fully po- 
larized electronic excitation, the simple protocol as described on 
the first page is applied (solid); lifetimes as below. (Dashed) Here 
the initial triplet populations based on our experimental findings 
(see text) are /?_ = 0.49, po = 0.02, p+ — 0.49 with associated 
lifetimes T_ = 0.57 ms, To = 0.02 ms, T + = 0.57 ms. In this case, 
the protocol is applied twice in succession, first for the population 
in 1 7+) and then for the population in |71). 



molecules with established techniques. 
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